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Bacteria possess a diverse range of mechanisms for inhibiting competitors, 
including bacteriocins, tailocins, type VI secretion systems and 
contact-dependent inhibition (CDI). Why bacteria have evolved such a wide 


array of weapon systems remains a mystery. Here we develop an agent-based 
modelto compare short-range weapons that require cell-cell contact, with 
long-range weapons that rely on diffusion. Our model predicts that contact 
weapons are useful when an attacking strain is outnumbered, facilitating 
invasion and establishment. By contrast, ranged weapons tend to be 
effective only when attackers are abundant. We test our predictions with the 
opportunistic pathogen Pseudomonas aeruginosa, which naturally carries 
multiple weapons, including CDI and diffusing tailocins. As predicted, 
short-range CDI can function at low and high frequencies, whilelong-range 
tailocins require high frequency and cell density to function effectively. 
Head-to-head competition experiments with the two weapon types further 
support our predictions: a tailocin attacker defeats CDI only when itis 
numerically dominant, but then we find it can be devastating. Finally, we 
show that the two weapons work well together when one strain employs 
both. We conclude that short- and long-range weapons serve different 
functions and allow bacteria to fight both as individuals and as a group. 


One of the most striking illustrations of Darwin's 'struggle for exist- 
ence" is the evolution of weaponry”’. Weapons-traits that evolved to 
injure and harm competitors—have evolved many times in animals, 
with examples in groups as diverseas trilobites, insects, mammals and 
dinosaurs’. Bacteria are a second group of organisms that commonly 
evolve weapons?^. Many clinical antibiotics were first isolated from 
bacteria that release them into the environment to inhibit competi- 
tors”. Ribosomally synthesized bacteriocins are deployed in a simi- 
lar manner and include both chemical toxins and phage-tail-derived 
tailocins, which physically punch holes in competitors". Bacteria also 
deploy close-range weapons that require contact between cells. Exam- 
ples include type VI secretion systems (T6SSs), which fire toxin-laden 
needles into competing cells", and contact-dependent inhibition 
(CDI) systems, which are toxin-loaded filaments anchored to the out- 
side of the cell”. The diversity of weapons seen in bacteria, therefore, 


certainly rivals that seen in animals. However, there is a notable dif- 
ference between the two groups. Excluding teeth and claws, whose 
primary evolutionary function is feeding, animals tend to carry a single 
weapon type: for example, horns or antlers or tusks”. Some animals 
are known to possess multiple weapons-the dinosaur Ankylosaurus 
magniventris had both horns and a bludgeoning tail club^*—but such 
examples appear to bethe exception. By contrast, bacteriacommonly 
carry multiple types of weapon’. 

Pseudomonas aeruginosa is a problematic opportunistic patho- 
gen, dueto its ability to withstand numerous antibiotics". Alongside 
its defensive capacity, this species is a striking illustration of how many 
weapons bacteria can carry. P. aeruginosa produces multiple bacteri- 
ocins and toxic small molecules, which serve as long-range weapons”. 
In addition, it can deploy CDI and up to three T6SSs as short-range 
weapons”. More generally, among species whose weapons have been 
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characterized in detail, many carry both short- and long-range weap- 
ons, including strains of Bacteroides fragilis'*"’, Pectobacterium caro- 
tovorum””, Burkholderia cepacia”, Chromobacterium violaceum™» 
and Myxococcus xanthus*°”’. What is the evolutionary basis for the 
prevalence of multiple weapons? One argument is that bacteria are 
simply more aggressive than other species such as animals, and that 
this favours the simultaneous use of multiple weapons. Consistent with 
this hypothesis, experiments suggest that bacteria engage in combat 
much more regularly than animals’, which typically avoid using their 
weapons’. However, a general increase in aggression does not explain 
why bacteria carry multiple types of weapon, as opposed to simply 
just investing more ina single type. We hypothesized that bacteria 
carry multiple weapons because they serve different functions during 
competition. We further reasoned that this explanation is most com- 
pelling for weapons that function at different ranges, a factor that can 
strongly influence the outcome of bacterial contests”. We therefore 
sought to test this hypothesis by performing a direct comparison of 
the competitive benefits of short- versus long-range weapons. 

Wefirstemploy arealistic agent-based modelof bacterial compe- 
tition that has previously been used to understand the evolutionary 
function of single weapons?"??, The power of this framework is that 
it allows one to rapidly study a wide variety of competition scenarios 
with relative ease, while being realistic enough to generate focused 
predictions for empirical testing. The model predicts that short- and 
long-range weapons do indeed have the potential to serve different 
functions. We test these predictions by genome editing P. aeruginosa 
strain PAO1 to generate strains that are susceptible to its own short- and 
long-range weapons (CDI and tailocins, respectively). This approach 
allows us to directly compare weapons' functioning in a controlled 
genetic background, and thereby investigate the relative benefits of 
thetwo weapontypes. Insupport ofthe modelling, we find that short- 
and long-range weapons can provide different advantages during 
combat. Contact weapons remain effective when an attacking strain 
is outnumbered, facilitating invasion and establishment. By contrast, 
ranged weapons are most effective when attackers are abundant, but 
herethey proveto be a devastating form of attack. 


Results 

An agent-based model of short- and long-range weapons 
Wefirst employ an established computational model where each cell 
is simulated as an individual agent (hence, ‘agent’ or ‘individual based’ 
model; Methods)? ?. Bacterial cells are seeded onto a two-dimensional 
(2D) surface where they grow and divide to fill a vertical space, as can 
occur in a biofilm for example, or at the intestinal mucosa of a mam- 
malian host. Cells interact with each other both physically—pushing 
and displacing each other as they grow and collide—and chemically, by 
producingtoxins that can inhibit strains of a different genotype (Fig. 1 
andSupplementary Video 1). When cells reach a specified height they 
are removed, mimicking dispersal or shedding from the top of the 
community. Bacterial weapons can vary substantially in a wide range 
of properties, including method of delivery (short- or long-range), 
quantity produced, deadliness, cost of production and, of course, 
thestrainor species where they are found". This can make like-for-like 
comparisons of different weapons empirically challenging. However, 
with a model one can precisely define, and systematically vary, such 
properties of bacterial weaponsto study general principles (Methods). 
The modelis also spatially explicit, which is particularly important for 
contact-based weapons whose action depends upon the occurrence 
of physical contact between cells??^?, 


Modelling predicts distinct strengths of short- and long-range 
weapons 

We use our model to compare short- and long-range bacterial weap- 
ons that—other than their range of effect—are as similar as possible. 
Attacking cells can 'fire' ashort-range weapon at random from their cell 


surface, intoxicating any susceptible cell(s) that are contacted (Fig. 1a), 
which is intended to simulate contact weapons including T6SSs???? and 
CDI”. Alternatively, attacking cells can release a diffusible factor into 
the environment, which could represent a range of toxins, including 
asmall-molecule antibiotic, a ribosomally synthesized bacteriocin 
or atailocin. The rate at which toxins are exported out of attacking 
cells is matched for short- and long-range weapons, and is controlled 
via a secretion rate parameter, Ke» Toxins also have matched poten- 
cies: both contact and diffusible toxins are lethal once intracellular 
concentrations exceed a set threshold T.. We assume that producing 
either toxin incurs an equivalent growth-rate cost in attacking cells 
that is proportional to the secretion rate: that is, increasing the toxin 
secretion rate results in a lower growth rate (Fig. 1b and Methods)^*^. 
Our goal is to compare weapons that differ only in their range, while 
keeping all other features identical. As a result, we assume that the cost 
of useis equivalent for both short- and long-range weapons. 

We begin by modelling competitions between an attacker strain, 
which either has a short- or a long-range weapon, and a second strain 
that is susceptible to the weapon. We look at a wide range of competi- 
tion scenarios, varying the initial frequency of the attacker, initial 
density of cells and the amount that the attacker invests in its weapon 
(that is, toxin secretion rate). In each case, the two strains are allowed to 
grow and interact for a set period of time (10 h; Fig. 1c), after which we 
compare the final attacker:susceptible cell ratio to its initial value (‘com- 
petitive advantage’; Methods). Both weapons tend to perform better 
when cells are seeded at high initial density, because this promotes 
cell-cell contact for short-range weapons” and toxin accumulation 
for long-range weapons. Nevertheless, there remains a clear difference 
between the two weapon types. Across the majority of scenarios tested, 
the contact-dependent weapon provided an advantage to the attacker 
(Fig. 1d and Extended Data Fig. 1a). By contrast, the long-range weapon 
is more sensitive to starting conditions, requiring a higher secretion 
rate to be effective (Supplementary Fig. 1b) and reaching its full poten- 
tial only when the attacker is seeded at high frequency and high density. 

The models suggest a particular advantage of contact-dependent 
weapons: they are effective even when users are at a numerical disad- 
vantage. This benefit is likely to be most substantial when invading an 
established population, because here invaders will be outnumbered 
by residents. We explored this scenario further by modelling estab- 
lished biofilms of susceptible cells and simulating the late arrival of 
attacker cells (by replacing some of the susceptible cells at random 
with attackers). In this scenario, cells using the short-range weapon 
were able to successfully invade an established population, where 
increasing toxin secretion rate increased their success (Fig. le,f and 
Extended Data Fig. 1c). Conversely, long-range weapons never enabled 
invasion, as producer cells could never amass in sufficient numbers to 
kill the susceptible strain. 

In summary, even though we closely matched the properties of 
the two weapon types, the models predict that they perform differ- 
ently across the variety of competition scenarios we tested. In general, 
contact weapons work well across a range of frequencies, including 
cases whenastrainis rare, allowing it to invade. By contrast, long-range 
weapons function well only when a producing strain is relatively 
abundant, but here they are an effective form of attack. 


Using genome editing to generate strains for weapon 
comparisons 

Our modelling predicts that a long-range weapon will perform poorly 
atlow attacker frequency, which is consistent with the findings of sev- 
eral previous studies—both theoretical and empirical -showing that 
toxin production is most effective when attackers are abundant??? *!, 
However, to test our predictions on the relative benefits of short- versus 
long-range weapons, a well-controlled comparison of the two types of 
weaponis required. To do this, we turned to the opportunistic pathogen 
P. aeruginosa (strain PAOT), which naturally carries both short- and 
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Fig. 1| Agent-based modelling predicts that contact weapons are more 
robust to changes in frequency, density and secretion rate. a, Contact toxins 
(top): producing cells can deliver toxins to neighbouring cells. If asusceptible 
cell (yellow) is within range, the toxin is injected (left dashed circle) and the 
susceptible cell dies; otherwise the toxin is wasted (right dashed circle). 
Diffusing toxins (bottom): when the local concentration of a diffusible toxin 
exceeds a threshold (within dashed line), susceptible cells die. b, Cells secrete 
toxins, incurring a growth-rate penalty proportional to the amount of toxin 
being secreted (secretion rate). c, Snapshots of competition outcomes for 
attackers with contact-dependent toxins (blue cells, left column) or diffusible 
toxins (magenta cells, right column). Unarmed susceptibles (yellow cells) die 
upon lethal toxin exposure (black cells). The contact weapon performs better 
atlower frequencies than the diffusible weapon. Snapshots show cropped 
(150 um) sections of the 300-um-wide, 2D simulation domain; below the black 
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line (arrow) represents the lethal concentration for the diffusible toxin. Scale bar, 
50 um; inoculum, 100 cells. d, Quantification of competition outcomes for two 
initial cells densities: ‘low’ (10 cells inoculum) and ‘high’ (100 cells inoculum). 
Competitive advantage assesses the fold change in the attacker strain compared 
with its competitor from the beginning to end of the simulation (Methods). 
Horizontal bars indicate the mean from multiple simulations (n = 6). e, Snapshots 
of competition outcomes for invasion scenario (invader frequency: 10%), with 
contact-dependent and diffusible-toxin-armed attackers coloured as in c. Scale 
bar, 10 um. Successive timepoints (rows) show fates of initially rare attackers 
following random inoculation into confluent biofilms of susceptible cells. 

f, Quantification of competition outcomes for invasions (invader frequency: 10%) 
using the same competitive advantage metric as in d, quantified as a function 

of toxin secretion rate. Horizontal bars indicate the mean from independent 
simulations (n = 6). 


long-range weapons. Of its long-range weapons, the mechanical tai- 
locins (specifically its R-type pyocins) are known to be highly effective 
weapons under biofilm-like conditions^"*. For short-range weapons, 
P. aeruginosa has a dedicated antibacterial T6SS", but it is under 
complex regulation and typically fires only in response to incoming 
attacks??^*^, Therefore we instead chose to focus on a CDI system of 
P.aeruginosa'* as a short-range weapon to test predictions. The fitness 
costs of these weapons potentially differ. Tailocins are induced by DNA 
damage and released through the self-lysis of a sub-population of cells”. 


CDI systems are regulated by unknown signals andare expected to be 
produced by all cells of a population, but their use does not require 
cell death*’. The regulatory strategies of weapons can have important 
effects on their effectiveness*®, which we do not study here. Never- 
theless, both weapons have been shown to provide clear advantages 
to their users, which made them a good choice for representative 
short- and long-range weapons. 

For both weapons, we used genome editing to generate a strain 
that is susceptible to the weapon but otherwise well-matched to the 
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attacker, allowing us to study the effects of each weapon on bacterial 
competition. For CDI, this was straightforward: deleting the three gene 
locus that encodes the CDI transporter, toxin filament and immunity 
(PAO040-PA0041), resulting in a strain that does not have CDI and is 
susceptible to CDI. For the tailocins, we selected pyocin R2, but here 
resistance is more complex as it is determined by the composition of 
thelipopolysaccharide (LPS) moieties of the outer membrane”. Here, 
we engineered a susceptible strain by deleting both the pyocin R2 
locus and the gene wbpL, which causes a deficiency in the LPS that in 
our strain background leads to susceptibility to pyocin R2 (Extended 
Data Fig. 2). However, this deficiency in the LPS caused a competitive 
disadvantage to the susceptible strain, independent of the effects of 
tailocins (Extended Data Fig. 3). To generate a near-matched attacker 
strain, therefore we made a AwapR deletion in the wild-type strain, 
which causes a similar LPS deficiency but not one that leads to pyocin 
R2 susceptibility (Extended Data Fig. 2 and Methods). 

In the absence of the weapon-mediated advantage of pyocin R2 
(hereafter ‘tailocin’), the wapR mutation puts the attacker at a mod- 
erate disadvantage relative to the wbpL mutation in the susceptible 
strain (Extended Data Fig. 4). We quantified this difference and used it 
to adjust the predicted competitive advantage of the attacker strain in 
our experiments, in order to estimate the effects of the weapon alone 
(Extended Data Fig. 4 and Methods). However, in practice this adjust- 
ment has little impact on the data because the benefits of the tailocin, 
when they are seen, massively outweigh this moderate cost of wapR 
deletion (for further discussion of these LPS mutations, see Methods). 


Experiments with P. aeruginosa show distinct advantages to 
short- and long-range weapons 

With these strains, we could then test our modelling predictions for 
a scenario where attackers armed with either a contact or diffusible 
weapon compete against a susceptible strain (Methods). Both weapons 
are probably influenced by the spatial structuring of the environment. 
With structured environments, for example, there is the potential for 
diffusion limitation that might constrain long-range weapons. Spatial 
structure can also limit short-range weapons if single genotypes are 
growing in distinct patches, because this can reduce contacts between 
attacker and target cells??^?. The spatial structure of the environment, 
therefore, may affect short- and long-range weapons differently. 
To include the potential for such effects, we performed competitions 
when cells are growing on agar (the ‘colony biofilm model)? ?, This 
assay allows us to capture the dense, spatially structured conditions 
thought to be typical of bacterial communities”, and both weapons are 
expected to function well in this context (Methods). The interior and 
edge of bacterial colonies represent distinct competition scenarios, due 
tothe much greater potential for population expansion at the edge”. 
We therefore decided to sample each region separately (Methods), 
although the competition outcomes show similar trends between the 
two. Consistent with previous work, we find that both weapon systems 
have the potential to provide large competitive benefits for an attacking 
strain*?“°>°°”, As in the models, outcomes were quantified by compar- 
ing the final ratio of attacker:susceptible cells with its initial value (the 
‘competitive advantage’). 

Overall, we find good support for our modelling predictions, 
despite the potential for differences in the cost of using tailocins and 
CDI. As the models predict, high attacker frequency is most important 
for the effectiveness of the long-range weapon (tailocin; Fig. 2b). In 
comparison, the short-range weapon (CDI) almost always performs 
equivalently or better at intermediate and low frequencies, especially 
at the colony edge (Fig. 2c). The CDI experiments also include cases 
where weapon performance peaks at intermediate frequency asinthe 
model (compare Fig. 1d, high density, with Fig. 2b, density 10* cells). 
For the tailocin competition, we observe an improvement in weapon 
performance as initial cell densities increase. This pattern is again 
predicted by the model (Fig. 1), although the effect is substantially 


stronger in the experiments, which cover much larger total numbers 
and ranges of initial cell density than are possible with the modelling. 
We note that previous work on long-range inhibition did not find this 
density-dependence?, but the inhibitory mechanism in this study 
relied on the diffusion of a quorum-sensing signal and a synthetic 
gene circuit, not a bacterial weapon, so itis not directly comparable. 

Fluorescence micrographs ofthe colonies show how the competi- 
tions play out in space (Fig. 2a and Extended Data Fig. 5), and allow a 
finer assessment of spatial structure than can be achieved with sam- 
pling of the edge and the middle for counting. These images reveal a 
pattern wherethetailocin performs worse atthe very edge ofthe colony 
as compared to CDI (Fig. 2a middle bottom). This pattern is consistent 
again with the requirements for tailocins to build up to be effective, as 
this build up is expected to happen first in the colony interior and last at 
the colony edge. The genotypic patchiness seen in these images is also 
likely to help explain the variability in the quantification data. Patchi- 
ness, especially at low density, is more prominent at the edge, leading 
to variation between replicates due to the coarse nature of the sampling 
procedure. Nevertheless, the overall patterns and difference between 
weapon types remain clear and statistically significant (Fig. 2b,c). 

In summary, our experiments show that both weapons can be 
highly effective, but that the two weapons perform best under different 
conditions. Tailocins are extremely effective at high frequencies and 
densities, while CDI performs more consistently across conditions, 
including the unique ability to provide a competitive advantage when 
astrain starts out rare and at low density. 


Head-to-head contests of short- and long-range weapons 
Our first experiments examine the performance of each weapon type 
against susceptible cells that do not fight back. We next explore the case 
where users of the two weapons meet. In this situation, a weapon can 
potentially take on new significance, as eliminating competitor cells 
alsoserves to reduce incoming attacks. The agent-based model again 
predicts that both initial frequency and cell density can be critical to 
the performance of the weapons (Fig. 3a,b, Extended Data Fig. 6 and 
Supplementary Video 2). At low cell density, the impact of both weapons 
is limited, but the contact weapon user does gain an advantage when 
it starts in the majority. At high cell density, the long-range weapon 
user can win but, as in the single weapon competitions (Fig. 1d), this 
requires it to start at high frequency. When the long-range strain is at 
low or equal initial frequency, the contact weapon performs the best. 
Totestthese predictions, we competed astrain susceptibleto CDI 
against a strain susceptible to tailocins (pyocin R2; Fig. 3c,d). The CDI 
attacker was thus an LPS mutant (AwbpL; see above and Methods), which 
renders it susceptible to the tailocin. As this mutation causes changes 
to the cell envelope, we first checked that the CDI system remained 
functional. The gene deletion (AwbpL) did not prevent CDI from func- 
tioning, but it did reduce the advantage provided by CDI at lower initial 
densities (Extended Data Fig. 7). Going forward, therefore, we only 
provide data from the two highest initial cell densities (10° or 10° cells 
ul). We also focus on data from the interior of the colony going forward, 
as they are most reflective of typical biofilm growth”, but data from 
the colony edge are similar (Extended Data Fig. 8). As predicted by the 
high-cell-density model, only the CDI-using strain gains a significant 
advantage in equal-frequency competitions (Fig. 3d). Moreover, again 
as predicted, both weapon users are able to win when they start in the 
majority. This effect is particularly strong for the tailocin, which pro- 
vides a large competitive advantage when starting from high frequency. 


Contact and diffusible weapons are complementary 

Our findings thus show that short- and long-range weapons can pro- 
vide distinct advantages, which helps to explain why bacteria would 
carry both types. However, it is possible that these advantages are not 
provided simultaneously in practice. Therefore, we sought to experi- 
mentally study a single strain using both weapons, compared to each 
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Fig. 2| Experiments show the importance of high density and high frequency One-way ANOVA showed that initial density and ratio significantly affected both 
for long-range weapons. Colony competitions with P. aeruginosa PAO1 between weapons at the edge (CDI, density: P = 6.54 x 10 *, n = 80; CDI, ratio: P= 2.29 x 


wild-type and mutants susceptible to either CDI (short range) or pyocin R2 10%; tailocin, density: P= 1.89 x 105, n= 92; tailocin, ratio: P=1.43 x 10°). For b 
(tailocin, long range) inoculated from different densities (mean inoculum density and c, competitions were assessed via counts of CFUs. Competitive advantage 

1.9 x 10°, 10*, 105, 10° CFU pI"). a, Representative microscopy images from equal- assesses the fold change in the attacker strain compared with its competitor from 
frequency (1:1) competitions after 48 h of growth. All strains express constitutive thebeginning to end ofthe competition. The tailocin attacker advantage inb 
fluorescent protein genes and are false-coloured either blue (CDI attacker, top), and chas been adjusted for a disadvantage in the background genotype ofthe 
magenta (tailocin attacker, bottom) or yellow (susceptible, top and bottom). attacking strain (see Methods and Extended Data Figs. 3 and 4). Horizontal bars 
Scale bar, 500 um. b, Quantification of competition outcomes at the colony indicate the mean of independent biological replicates (n > 6; see Supplementary 
centre. One-way ANOVA showed that initial density and ratio significantly Table 5for exact values of n). Top brackets indicate a significant difference 
affected both weapons in the centre (CDI, density: P=1.45 x 10 5, n= 80; CDI, between the weapons (two-sided Welch's t-test, P < 0.05, Benjamini-Hochberg 


ratio: P=1.45 x 10 5; tailocin, density: P = 3.68 x 10 5, n = 91; tailocin, ratio: P = 5.01 correction for multiple testing; see Supplementary Table 5 for exact P values). 
x10). c, Quantification of colony competition outcomes at the colony edge. 
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genotype of the tailocin-using, CDI-susceptible strain (magenta) is AwapRACDI. 


weapon being used alone. To allow susceptibility to tailocins and make 
the comparison as fair as possible, we put all attacker strains in the 
same LPS genetic background (AwapR), and all susceptible strains in 
the AwbpL background. All competition outcomes are again adjusted 
forthe moderate disadvantage caused by the AwapR mutation, and we 
focus on high initial cell densities as before (see above and Extended 
Data Fig. 9). As predicted by the study of each weapon individually, 
across the competitions and contexts, the two weapons function in a 
complementary fashion, providing an equivalent or greater benefit 
than either weapon alone (Fig. 4 and Extended Data Figs. 9 and 10). 
That is, we find that using both CDI and tailocins can allow a strain to 
receive benefits from both. 


Discussion 

Bacteria use a vast variety of weapons to inhibit and kill competitors. 
Here we have shown that two major categories of weapon-short- and 
long-range—can provide distinct and complementary advantages 
to bacteria. We find that weapons that rely on contact between cells 
are generally effective whether a strain is rare or common, which is 


consistent with previous work on the CDI systems of Escherichia coli”. 
Conversely, weapons that diffuse across long ranges are more reliant 
on high producer cell frequency and density. However, under these 
conditions, long-range weapons can be particularly powerful in elimi- 
natingboth armed and unarmed competitors. These observations help 
to answer two related questions. First, why do short- and long-range 
weapons exist at all in bacteria? Here, our work suggests that a strain 
may benefit more from one weapon type or another, depending on its 
ecology. For example, if bacterial fitness is more determined by the 
ability of a given strain to invade communities than persist in a com- 
munity, contact-dependent weapons may be most useful. This matches 
with observations showing that pathogens including Salmonella? and 
Vibrio cholera? use the T6SS during invasion of established communi- 
ties, although gut-resident Bacteroides fragilis can also useits T6SS to 
repel invaders” and the HI-T6SS of P. aeruginosa fires only in response 
toincoming attacks", underlining the flexibility of short-range weap- 
ons (Figs. 1-3). 

The second question is this: why do many bacteria carry both 
short- andlong-range weapons? Again, the fact that the two weapon 
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Fig. 4 | The benefits of short and long-range weapons combine 

positively in P. aeruginosa. Quantification of competition outcomes in 

the colony centre for two initial cell densities (mean inoculum density 

1.9 x 10°, 10° CFU pI"). Competitive advantage assesses the fold change 

in the attacker strain compared with its competitor from the beginning 

toendofthe competition. Competitions where the attacker has just 

CDI (blue, left), just tailocins (magenta, centre) or both weapons 

(purple, right) show the advantage gained from using two weapons together 

as compared to just one. Data are adjusted to account for differences in 

competitiveness of the strain backgrounds (AwapR relative to AwbpL; see 

Methods and Extended Data Fig. 2). Horizontal bars indicate the mean from 

independent biological replicates, n 2 6; see Supplementary Table 5 for 

exact values of n). Stars above the double weapon data indicate a significant 

difference between the combination of weapons and either single weapon 

(blue and magenta), or just CDI (blue) (two-sided Welch's t-test, P « 0.05, 

Benjamini-Hochberg correction for multiple testing; see Supplementary 

Table 5 for exact P values). 


typescan perform distinct functions provides an answer to why a cell 
would carry both. However, the possession of multiple weapons may 
also carry other non-mutually exclusive advantages. If two weapons 
are simply more potent than one, this can also favour carrying several 
of them”. In addition, if a given competitor is likely to carry resist- 
ancetosome ofthe available weapons, natural selection may favour 
using multiple weapons simply to increase the likelihood that the 
competitor is susceptible to at least one. This explanation may be 
most important for cases where bacteria carry several of the same 
weapon type and functional differences are less pronounced, for 
example, P. aeruginosa releasing multiple S-type pyocin protein toxins 
and tailocins simultaneously via cell lysis'é. If the use of short- and 
long-range weapons can be so beneficial, why is it not more widely 
seeninorganisms otherthan bacteria? Short-range weapons, like CDI 
and T6SSs, can be used against competitors by lone cells and remain 
effective at low densities. However, the effectiveness of long-range 
weapons rests upon the strength in numbers that comes with group 
living. Mounting evidence suggests that bacteria commonly live, 
and fight, as part of large groups’. Our work suggests that it is this 
propensity for group living that has led to the widespread evolution 
of both short- and long-range weapons. 


Methods 

Agent-based modelling 

A major challenge in comparing bacterial weapons is that even within 
the broad categories of contact and diffusible, their method of deploy- 
ment varies substantially. The T6SS can be fired at random intervals 
to deliver toxins or only in response to incoming attack**, whereas 
CDI filaments are produced to decorate attacking cells in unknown 
numbers, but toxin translocation occurs only in responseto receptor 
binding on a target cell. The small-molecule antibiotics of Strepto- 
myces are secreted from intact cells, whereas F. coli's colicins and the 
pyocins of P. aeruginosa require cell lysis for release. These examples 
also demonstratethe high variability inthe cost ofusing weapons. The 
useof computational models is very amenableto these challenges, as 
we can unify the production of both contact and diffusible weapons 
under asingle parameter 'secretion rate' and maketheir growth costs 
(per unit secretion) equal. 

The behaviours of cellular groups are often emergent, in the sense 
that they can only be understood interms of collective effects that arise 
whenever there are many interacting organisms”. A strength of our 
modelling approach is that it is able to capture emergent and spatial 
effects, such as the importance of cell shape for bacterial competi- 
tion??, the importance of lytic toxins for T6SS attack”? and the need 
for strong reciprocation when the T6SS is used in response to incoming 
attacks”. For each of these examples, the models predicted novel biol- 
ogy that was subsequently validated empirically using experimental 
work????^5 giving us confidence in our modelling approach. 

All simulations were carried out using CellModeller?, an 
open-source software platform for running agent-based models of bac- 
terial growth. To model contact-based and diffusible toxin secretion, 
weimplement additional python-based CellModeller modules, whose 
source codes are available at https://github.com/WilliamPJSmith/ 
CellModeller. The key processes incorporated in our model aresumma- 
rized below; model variables and parameters are summarized respec- 
tively in Supplementary Tables 1 and 2. 


Model description. Cell growth and division. In our simulations, 
bacterial cells are represented as short capsules with a fixed radius 
R=0.5 um, and a birth segment length of 0.01 pm (equivalent to a birth 
volume V, = 0.54 jum?). Cells grow through elongation, dividing after 
doubling their birth volume Vo, plus a small random noise term €, = U(O, 
0.05) Vp. Each simulation timestep At, each cell grows by an amount 
proportional to that cell's current volume, V = KV, discretized as 
Vit + At) = V(t) (1+ K At). Here, K,.,, (s) represents the net per capita 
growth rate. The production of ranged or contact toxins is assumed to 
be costly, such that weapon users suffer a growth penalty proportional 
to their toxin secretion rate, Kgrow = Kmax(1-CKsec), With c, kmax and Ksec the 
prorata toxin cost, maximum growth rate and equivalent contact toxin 
secretion rate, respectively. For simplicity and computational expedi- 
ency, we assume throughout that nutrient access is non-limiting, such 
that kmax is independent of cells’ positioning ina community. 


Mechanical interactions. Mechanically, cells are modelled as rigid, 
elastic particles that push on one another as they grow and divide. 
Each simulation timestep, immediately following the growth stage 
outlined above, an energy penalty method is used to compute cell 
movements necessary to minimize total cell-cell overlap, subject to 
viscous drag forces acting on each cell. This process, described previ- 
ously in detail?^*5*, approximates the elastic repulsion forces acting 
between cells in physical contact. 


Contact-dependent toxins. As in previous publications????, we 
use a custom Python module to represent cell-cell antagonism via 
contact-dependent toxins. While this module was previously used to 
study T6SS-mediated interference, it is a generic representation of 
contact-dependent warfare that allows usto make general predictions 
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that should apply to many mechanisms, including CDI. Each simulation 
timestep, cells armed with contact toxins fire needles of length R, pro- 
jecting orthogonally from randomly chosen sites on their cell surface. 
The number of secretion events per cell per unit time is drawn froma 
Poisson distribution, whose mean is the secretion rate Kec After firing, 
each needle is checked to determine if it comes into contact with any 
other cell in the population (line-segment method). Successful hits 
are logged for each target cell and result in cell death if their number 
(excluding hits by kin cells) exceeds a lethal threshold N,;, = 1 (ref. 65). 


Diffusible toxins. We assume toxins to be freely diffusible solutes 
that kill susceptible cells when their local concentration u+ exceeds a 
lethal threshold T, (controlled by N,its; see below). To represent natural 
variability in toxin susceptibility“, lethal toxin threshold is drawn for 
each cell from a normal distribution, N(1, 0.2) at birth (we ignore this 
stochasticity in the discrete contact toxin model, because with mean 
Nits = L, the chances of any cell surviving more than one hit are approxi- 
mately 1:130,000). To model the toxin concentration field u, = uj(x, y) 
(kg; m?) for a given cell configuration, we use the reaction-diffusion 
equation 0u,/dt = D; V?u, + k,ap¢(x, y). Here, Dr, kr, a, p and $(x, y) are 
respectively the toxin diffusivity (m ?), the specific toxin production 
rate (s), the toxin yield per unit cell biomass (kg; kg; !], the cell biomass 
density (kg, m?) and the cell volume fraction function (unitless). In 
non-dimensional form, pseudo-steady-state solutions to this equation 
are given by V?u; = D; (x, y). The behaviour of this equation is governed 
by asingle parameter grouping, the Damkóhler number D; = Pkyap/D;T.. 
Conceptually, we vary D, by changing the toxin production parameter 
k,, such that increases in production always incur proportional increases 
in production cost. We compute pseudo-steady-state solutions with the 
finite element method, using the FEniCS Python library and supporting 
CellModeller modules. Solutions are evaluated on a 2D rectangular 
domain of dimensions L, by L, with crossed mesh element size A = 5 uum. 
Solutions are subject to mixed boundary conditions: periodic boundary 
conditions (left and right edges), Neumann boundary conditions (base 
edge) and Dirichlet boundary conditions (u, = 0 along top edge). Asin 
previous studies”, we assume that toxin is not subject to degradation 
or removal as part of its activity; toxin is only lost from the domain via 
leakage along the top edge. 


Parity between diffusible and contact-dependent toxins. Our model 
aims to compare diffusible and contact weapons ina like-for-like man- 
ner. For this comparison to be as fair as possible, we assume in all cases 
that both weapons involve the secretion of the same (hypothetical) 
toxin, at the same rate kyen with the same potency (lethal concentra- 
tion) Te, and at the same growth cost c. For simplicity, the secretion 
rate is fixed within each simulation, while in reality many bacteria 
responsively upregulate toxin production via competition sensing?! 
and other mechanisms”. Here we outline the equations that link these 
properties between the two models, bridging the (discrete) contact 
model with the (continuum) diffusible toxin model. For contact 
weapons, the per-cell secretion rate kz... is given by Kr = KeecT sec with 
Ke (nh) the per-cell secretion rate, and T.se (kgr) the mass of toxin 
released by each secretion event. For diffusible toxins, this can be 
expressed as ky cen = Kyapf? (terms defined as above) with P approximat- 
ing the cell volume. The contact-dependent secretion rate K,.. can there- 
fore be related to the diffusible toxin Damkóhler number, D,, as D; = 
T. KD, Te. To relate the potencies of contact and diffusible toxins, 
we assume that NM, is equivalent to the minimum number of contact 
events required to raise the intracellular toxin concentration to the 
lethal threshold Tc: My, = TcP/T;.., With Papproximatingacell's volume 
as before. Moreover, by combining these equations, we can eliminate 
the (unknown) contact toxin load 7... and compute D, = kee? /N pits) 
for the two weapons operating at equivalent secretion rates. This 
relation highlights that toxin diffusivity D; is a crucial parameter for 
our model, because increasing D, is equivalent to reducing effective 


diffusible toxin production while keeping contact toxin production the 
same. While previous work has explored the interplay between toxin 
diffusivity, cost and environmental structure, here we assume a fixed 
value for D, (that for colicin la”), whichis towards the lower limit of the 
range explored by previous work”. 


Simulation domain and protocols. All model simulations are run on 
an L,-by-L, rectangular domain with lateral periodic boundary condi- 
tions, representing a vertical 2D slice through a bacterial community. 
The base of the domain, y = 0, is an impenetrable substrate. The domain 
width L, is fixed; its height L, is set to track the maximum height p, of 
bacterial cell groups as L, = max(p,) + 6, with 6 the diffusive boundary 
layer thickness. To represent cell detachment through mechanical 
sloughing, cells are removed from the simulation once they reacha 
height Asiougner: This provides a simple representation ofa biofilm with 
alimited carrying capacity, and (by capping the number of simulated 
cells) allows simulations to be run for longer than if the population 
were allowed to increase indefinitely. All simulations are initiated by 
randomly scattering cells along the base edge of the simulation domain. 
Unless otherwise indicated, simulations run for a fixed duration of 
400 steps (10 h of simulated time, timestep At = 0.025 h). Invasion 
simulations involve a two-step process: first, the domain is inoculated 
as above using only susceptible cells, which are then allowed to divide 
until the cell group reaches the slougher height Agoughers Then, a random 
sub-population of the susceptible cells is converted to the invading cell 
type, and the simulation is allowed to run for up to 1,000 steps (25h), 
terminating early if either cell type is lost from the simulation. 


Computation and postprocessing. All agent-based simulations were 
run using a 2017 Apple MacBook Pro laptop computer, with concurrent 
simulations distributed between an Intel 3.1 GHz quadcore i7-7920HQ 
CPU, an Intel HD 630 Graphics card and an AMD Radeon Pro 560 Com- 
pute Engine. Simulation data were visualized using Paraview software 
(5.4.0), and analysed using custom MATLAB and R scripts. 


Experiments 

Strain construction. Deletion mutants were constructed using stand- 
ard two-step allelic exchange methods using the vector pEXG2 and 
Gibson assembly”. Strains constructed and used are summarized in 
Supplementary Table3. Primer sequences for up/downstream regions 
and exterior confirmation primers are listed in Supplementary Table 4. 
Constructed deletion vectors were introduced into P. aeruginosa PAO1 
by conjugation with £. coli JKE201 (ref. 76) and transconjugants were 
selected on lysogeny broth (LB) agar with 50 pg mI! gentamycin. After 
counter-selection on LB no salt, 10% sucrose, colony PCR positive clones 
were confirmed by Sanger sequencing (Source Bioscience) and genta- 
mycin sensitivity confirmed. Stocks in LB 2096 glycerol were stored at 
-80 °C. For strains with multiple deletions, they are listed in the order 
the deletions were made. Strains were subsequently constitutively 
tagged with eYFP and mScarlet (Sujatha Subramoni, unpublished) using 
pUCI8-mini-Tn7-GmR" delivered by conjugation with £. coli S17À and 
selected on Pseudomonas Isolation Agar with 100 ug ml! gentamycin. 


Engineering strains for weapon comparisons. For CDI-mediated 
competition, we deleted the three gene locus that encodes the CDI 
transporter, toxin filament and immunity (PAOO40-PA0041), which 
resulted ina strain that does not have CDI and is susceptible to CDI. For 
pyocin R2, we found that deletion of wbpL led to susceptibility. wbpL 
isatransferasethatinitiates the formation of LPS chains by adding the 
first sugar to the undecaprenol-phosphate carrier, and in a previous 
study a cosmid library-derived insertional inactivation mutant was 
found to be resistant to pyocin R2 (ref. 49). For our purposes, we made 
anin-frame deletion mutant of wbpL from scratch, which grew poorly 
inliquid media. We reasoned that this was dueto inhibition by its own 
pyocins. Consistent with this, making the wbpL mutant in a pyocin R2 
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deletion background restored growth. Moreover, this strain was then 
found to be susceptible to pyocin R2 from its parent. Finally, to confirm 
that wbpL was responsible for this susceptibility we complemented the 
in-frame deletion mutant of wbpL with a copy of wbpL on a plasmid, 
which restored pyocin resistance (Extended DataFig. 4). Onthis basis, 
we are confident that, in our strain background, deleting wbpL leads 
to pyocin susceptibility. 

wbpL was PCR-amplified from P. aeruginosa PAO1 genomic DNA 
and inserted into the expression vector pSEVA524 (Rubén de Dios, 
Eduardo Santero and Francisca Reyes-Ramírez, unpublished)" by Gib- 
sonassembly (NEB Hifi Assembly, NEB Location). After conjugation with 
E.coli JKE201andselection on LB10 ug ml” tetracycline a positive clone 
andaclonecarrying the empty vector were tested for susceptibility to 
pyocin R2. Strains were grown overnight at 37 °C in LB (with 10 pg mI 
tetracycline), then 1 ml mixed with 7 ml 0.75% LB agar and poured 
onto a warm LB plate to generate an overlay. Attacker cultures were 
also grown similarly and R pyocins were prepared by filter sterilizing 
supernatant from overnight cultures. Twenty microlitres was spotted 
onthe overlays, which were then dried and incubated overnight at 37 °C 
prior to photographing using a Gel Doc imager. 

We also made a deletion of wapR, which attaches the first 
L-rhamnose to the LPS core, initiating LPS capping, which did retain 
resistance to pyocin R2. We then had two strains, AwbpL and AwapR, 
that are LPS defective, but oneis susceptible to pyocin R2 and one is not. 
These strains were thus used for the long-range diffusible weapon exper- 
iments. To account for all fitness differences due to these LPS biosynthe- 
sis mutations inthe absence of effects from pyocin R2, AR2AwapR was 
competed against AR2ZAwbpL (Extended Data Fig. 4), and this difference 
was used as a baseline for comparison when the pyocin was present in 
the AwapR strain. All strains are available upon request. 


Culturing and the colony biofilm model. Strains were recovered from 
cryo stocks by streaking on LB 1.5% agar and incubating overnight at 
30 °C. LB 1.5% agar for competitions was prepared immediately prior 
to competition setup by pouring 20 ml into a Petri dish and allowing 
it to set for 15 min in a laminar flow hood. Colony competitions were 
prepared as previously by scraping cells off the overnight plate and 
resuspending cells to an initial OD600 of1 (ref. 51). Strains were mixed 
at defined ratios of 1:10, 1:1 and 10:1 then serially diluted 10-fold and 
1glspottedonthe prepared plate to generate competitions at various 
initial ratios and densities. Initial culture density was determined by 
serially diluting and spot plating. Independent biological replicates 
were performed with each attacker/susceptible combination carrying 
opposite fluorescent markers. 

We performed competitions when cells were growing on agar 
(the ‘colony biofilm model")? In nature, most bacteria live in densely 
packed communities, such as surface-associated biofilms, which are 
densely packed and spatially structured ^*^^, These high-density con- 
ditions are where both short- and long-range weapons are expected to 
function at their best, because they ensure plentiful cell-cell contacts 
and the potential for factors released in the cells to build up to high 
concentrations‘®*!*°*!, Bacterial weapons can also strongly influence 
the spatial structure of competing strains, and vice versa; effects that 
are not captured in liquid culture??9*??^* Finally, P. aeruginosa's CDI 
system is upregulated in structured static cultures, as compared to 
shaking culture, again suggesting these are the conditions where it 
has most impact^*. Growing bacteria on agar captures these dense 
and structured conditions in a highly tractable manner, allowing large 
numbers of competitions and conditions to be studied and imaged ^^ ^* 
(Fig. 2). Colonies also represent a good match to our model framework, 
which is again focused on high-density, biofilm-like conditions. 


Imaging of colonies and quantification of competition outcomes. 
After 48 h of growth at room temperature, colonies were imaged 
using a Zeiss Axio Zoom V16 microscope with a Zeiss MRm camera, 


0.5X PlanApo Z air objective and HXP 200 C fluorescence light source. 
Forty-eight hours allowed colonies inoculated from the lowest initial 
densities to grow sufficiently to be observed. All colonies from a sin- 
gle set of frequencies and densities were imaged at the same zoom 
(between x1and x2.5). To make the composite images shown in the 
figures, the display histograms of each channel were scaled to the 
minimum and maximum values found in the entire set of frequency and 
density, meaning images can be compared within sets (that is, weapon 
competitions) but not between. For Extended Data Fig. 5, all images 
were treated as a set, so comparisons can be made across all images. 
After microscopy imaging, colonies were sampled with a 10 pu pipette 
at boththe centre and edge ofthe colony into 0.9% saline. Samples were 
homogenized, serially diluted and 5 ul spotted onto LB or LB 50 ug mI! 
gentamycin and incubated at 30 °C overnight. Colonies were counted 
to determine the final ratio ofthe two strains. 


Calculation of competitive advantage. Using the initial density 
counted from the original inoculum cultures and the known inoculum 
ratios, the initial ratio of attacker:susceptible strains was determined. 
The final ratio was determined from the colony forming unit (CFU) 
counts of the serially diluted centre and edge samples, witha detection 
limit of 2,000 CFU mI" used to replace zeros and prevent dividing by 
zero. Competitive advantage was defined asthe final ratio/initial ratio 
andis plotted on log axes. 


Statistical tests. Data were checked for normality by inspecting 
histograms and quantile-quantile plots, then confirmed using the 
Shapiro-Wilk test. As most (27596 in each set) groups were normal, 
two-sided t-tests (Welch's) were used to find significant differences 
between weapons. To account for multiple hypothesis testing, each 
group oftests was corrected usingthe Benjamini-Hochberg method. 
Data analysis and figures were carried out in R version 3.6.3 using the 
following packages: tidyverse version 1.3.2 (ref. 85), dplyr version 1.0.9, 
broom version 1.0.1, ggplot2 version 3.3.6, ggfx version 1.0.1, patchwork 
version 1.1.2 and scales version1.2.1. Microscopy images were prepared 
for presentation using Image) version 1.530%. 


Reporting summary 
Further information on research design is available in the Nature Port- 
folio Reporting Summary linked to this article. 


Data availability 
All data are available at https://figshare.com/articles/dataset/Booth 
SmithFoster2023/23177156. Source data are provided with this paper. 


Code availability 

All code for the simulations was performed in a custom version of 
CellModeller 4.2.1, available at https://github.com/WilliamPJSmith/ 
CellModeller. Code for the data analyses is available with the data. 
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Extended Data Fig. 1| Agent-based modelling shows differences between 
weapons due to initial density of competitions, weapons depend differently 
on toxin secretion rate and that contact weapons better facilitate invasion 
than diffusible weapons at equivalent secretion rates. a Quantification of 
competition outcomes for all tested densities (secretion rate: 100). Densities of 
10 and 100 cells correspond respectively to ‘Low’ and ‘High’ starting densities 
showninFig.1.Competitive advantage assesses the fold change in the attacker 
strain compared to its competitor from the beginning to end ofthe simulation 
(Methods). Horizontal lines indicate the mean from multiple simulations 

(n = 6). b Outcomes of competition simulations over a range of secretion rates 


andinitial attacker frequencies (1096, 3096, 50?6, 7096, 90?6); initial density:150 
cells. Competitive advantage assesses the fold change in the attacker strain 
compared to its competitor from the beginning to end of the simulation 
(Methods). Horizontal lines indicate the mean from multiple simulations (n - 7). 
cQuantification of competition outcomes for invasions as a function of secretion 
rate. Invasion outcome is the same as competitive advantage (the fold change 
intheattacker strain compared to its competitor from the time of invasion to 

the end of the simulation). Horizontal lines indicate the mean from multiple 
simulations (n > 6, see Supplementary Table 5 for more details; only invasions 
wherethe invader was still present at the end ofthe simulation were analyzed). 
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Extended Data Fig. 2 | Lipopolysaccharide biosynthesis genes affect 
susceptibility to pyocin R2. Images of agar overlay assays showing pyocin R2 
zones of clearing (arrows) for different LPS mutants. The strain in the overlay 
isindicated in the center of each plate. The source strain for the pyocin R2 is 
indicated at the corners. Pyocins were prepared by sterile filtering supernatant 
fromovernight cultures. Overlays were prepared by mixing 1 mL of overnight 
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culture with 7 mL 0.75% LB agar then thoroughly drying. a AwapR shows no zones 
of clearing. b AR2 AwbpL (the entire pyocin R2 gene cassette was first deleted 
from this strain, then wbpL deleted second) shows zones of clearing from WT and 
AwapR, but not when pyocin R2 is deleted. c Complementing AR2ZAwbpL with 
empty vector pSEVA-524 does not rescue clearing. d Complementing AR2AwbpL 
with pSEVA-524 carrying wbpL shows no clearing from WT or AwapR. 
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Extended Data Fig. 3 | Lipopolysaccharide biosynthesis gene deletions affect areclosely matched and look similar to wild-type competed against AR2 (top). 
competition outcomes in the absence of pyocin R2. Microscopy showsthatthe Competitions were inoculated with -2 x 10° cells/uL. Images are representative 
AwbpL strain cannot compete with wild-type P. aeruginosa, evenintheabsenceof from three independent experiments. 

killing by tailocins (pyocin R2). Conversely, AwbpL and AwapR (second from top) 
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Extended DataFig. 4 | Deletion oflipopolysaccharide biosynthesis gene to end ofthe competition. Horizontal lines indicate the mean from biological 
wapR causes a disadvantage compared to deletion of wbpL when bothstrains replicates (n > 6, See Supplementary Table 5 for exact n values). The mean 

have pyocin R2 deleted. Quantification of colony competition outcomes (70.637) across all replicates from all densities and inoculum ratios for the center 
between lipopolysaccharide biosynthesis mutants in the absence of pyocin R2. was significantly different from O (One sided Welch's t-test, t = 20.48, 

Colonies were inoculated at the stated initial ratios and densities (meaninoculum df=111, p = 2.2e-16), so was used as the baseline advantage of AwbpL over AwapR. 
density 1.8 * 105, 10*, 10°, 10° CFU/LL). Competitive advantage assesses the fold This difference in advantage was subtracted from all competitions involving 
change in the attacker strain compared to its competitor from the beginning strains with these LPS biosynthesis gene deletions. 
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Extended Data Fig. 5| Colony competitions of contact dependent inhibition 
(CDI) and tailocins highlight differences between contact and diffusible 
toxins. Representative microscopy images of colony competitions inoculated 
from different starting densities (mean inoculum density 1.9 * 10°, 10*, 10°, 

10° CFU/LL). and initial ratios of attacker to susceptible cells taken after 48 h 

of growth. All strains are expressing constitutive fluorescent protein genes 

and false-coloured either blue (CDI attacker, top), magenta (tailocin attacker, 


10:1 


Initial Density 10° 


Initial Density 10° 


104 


Initial Densi 


Initial Density 10° 


500 um 


500 um 


1:10 


1:1 
Initial Attacker:Susceptible 


10:1 


© Susceptible Q CDI @ Tailocin 


bottom) or yellow (susceptible, top and bottom). Scale bar indicates 500 um. For 
the CDI competitions the attacker was wild-type and the susceptible has the CDI 
toxin and anti-toxin deleted. For the tailocin competitions, the attacker is AwapR 
and the susceptible strain is ARZAwbpL. Images shown here are representative, 
images were taken for every colony sampled (data presented in Fig. 2), exact 
values ofn are detailed in Supplementary Table 5. 
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in the simulation. Competitive advantage assesses the fold change in the attacker 
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Extended Data Fig. 7 | Outcomes of colony competitions shows that CDI 
remains functional in LPS biosynthesis gene mutants, but its effectiveness is 
diminished at low densities. Outcomes of CDI mediated competitions in wild- 
type (WT, green) or asymmetric LPS backgrounds. For these cases, both strains 
also havetailocins (pyocin R2) deleted. The attacking CDI+ strain is either AwapR 
against a CDI susceptible AwbpL (mauve) or CDI+ AwbpL against CDI susceptible 
AwapR (orange). Colonies were inoculated at the denoted initial densities (mean 
inoculum density 2.0 * 10?, 10*, 10°, 10° CFU/pL) and quantified by sampling, 


AwapR CDI-* vs AwbpL CDI- E AwbpL CDI+ vs AwapR CDI- 


plating and counting colony forming units after 48 h of growth. Horizontal lines 
indicate the mean from biological replicates (n > 4, see Supplementary Table 5 for 
exactn values). Competitive advantage assesses the fold change in the attacker 
strain compared to its competitor from the beginning to end of the competition. 
Top brackets indicate a significant difference between each single weapon and 
the combination of weapons (two-sided Welch's t-test, p « 0.05, Benjamini- 
Hochberg correction for multiple testing, see Supplementary Table 5 for 

exact p values). 
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Extended DataFig. 8 | Head-to-head weapon competitions between short Welch's t-test, p « 0.05, Benjamini-Hochberg correction for multiple testing, 
andlong-range weapon users at the colony edge. Quantification of direct see Supplementary Table 5 for exact p values). Stars indicate a competitive 
weapon colony competition outcomes at the colony edge by sampling, plating advantage significantly different from O (one-sided Welch's t-test, p « 0.05, 

and counting colony forming units. Competitive advantage assesses the fold Benjamini-Hochberg correction for multiple testing, see Supplementary Table 5 
change in the attacker strain compared to its competitor from the beginning to for exact p values). Competitions were inoculated with different initial densities 
end of the competition. Values above 1 (10°, dashed line) indicate an advantage (mean inoculum density 2.3 *10°, 10° CFU/LL). The genotype of the CDI using, 
for CDI, while values below 1 (that is 10 7, 10“) indicate an advantage for tailocins. tailocin susceptible strain (blue) is AR2ZAwbpL. The genotype ofthe tailocin 
Horizontal lines indicate the mean from biological replicates (n = 8). Top using, CDI susceptible strain is AwapRACDI. 
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Competitions wherethe attacker has just CDI (blue, left), just tailocins (magenta, Welch’s t-test, p « 0.05, Benjamini-Hochberg correction for multiple testing, see 
centre) or both weapons (purple, right) show the advantage gained from using Supplementary Table 5 for exact p values). Data from colony edge are noisier 
two weapons together as compared to just one. Data are adjusted to account thaninthe colony center but patterns are consistent with the colony interior. 
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Extended Data Fig. 10 | Microscopy images of colony competitions with fluorescent protein genes and false-coloured either purple (attacker) or yellow 
doubly-armed attackers. Representative microscopy images (taken after 48 h (CDI and tailocin susceptible). Scale bar indicates 500 um. The genotype ofthe 
ofgrowth) of colony competitions inoculated from different starting densities attacker is AwapR. The genotype of the susceptible strain is ACDIAR2ZAwbpL. 
(mean inoculum density 1.9 * 10°, 10° CFU/pL). and initial ratios of attacking Images shown hereare representative, images were taken for every colony 

and dual CDI/tailocin susceptible cells. All strains are expressing constitutive sampled (data presented in Fig. 4). 
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No statistical method was used to predetermine sample size. For all experiments, a minimum of six biological replicates were performed to 
allow for three replicates using each combination of fluorescent protein reporter tags. A minimum of 6 independent biological replicates for 
each condition was selected as a sample size as this was expected to cover the expected variability and the observation of clear trends shows 
this was true. 


No experimental data were excluded from analysis. 


Experiments were performed on separate days with independent pre-cultures, and replicated at least three times with both possible 
combinations of fluorescent protein reporter tags, resulting in clear trends indicating reproducibility. 


Allocation of samples was not random all cultures were explicitly set up using different combinations of strains to test the desired outcomes 
of competitions. 


Blinding is not relevant to this study as all experiments were set up under explicit conditions to test the desired outcomes of competitions 
between specific strains. 
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Briefly describe the study type including whether data are quantitative, qualitative, or mixed-methods (e.g. qualitative cross-sectional, 
quantitative experimental, mixed-methods case study). 


State the research sample (e.g. Harvard university undergraduates, villagers in rural India) and provide relevant demographic 
information (e.g. age, sex) and indicate whether the sample is representative. Provide a rationale for the study sample chosen. For 
studies involving existing datasets, please describe the dataset and source. 


Describe the sampling procedure (e.g. random, snowball, stratified, convenience). Describe the statistical methods that were used to 
predetermine sample size OR if no sample-size calculation was performed, describe how sample sizes were chosen and provide a 
rationale for why these sample sizes are sufficient. For qualitative data, please indicate whether data saturation was considered, and 
what criteria were used to decide that no further sampling was needed. 


Provide details about the data collection procedure, including the instruments or devices used to record the data (e.g. pen and paper, 
computer, eye tracker, video or audio equipment) whether anyone was present besides the participant(s) and the researcher, and 
whether the researcher was blind to experimental condition and/or the study hypothesis during data collection. 


Indicate the start and stop dates of data collection. If there is a gap between collection periods, state the dates for each sample 
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Data exclusions If no data were excluded from the analyses, state so OR if data were excluded, provide the exact number of exclusions and the 
rationale behind them, indicating whether exclusion criteria were pre-established. 


Non-participation State how many participants dropped out/declined participation and the reason(s) given OR provide response rate OR state that no 
participants dropped out/declined participation. 


Randomization If participants were not allocated into experimental groups, state so OR describe how participants were allocated to groups, and if 
allocation was not random, describe how covariates were controlled. 


Ecological, evolutionary & environmental sciences study design 


All studies must disclose on these points even when the disclosure is negative. 


Study description Briefly describe the study. For quantitative data include treatment factors and interactions, design structure (e.g. factorial, nested, 
hierarchical), nature and number of experimental units and replicates. 
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Research sample Describe the research sample (e.g. a group of tagged Passer domesticus, all Stenocereus thurberi within Organ Pipe Cactus National 
Monument), and provide a rationale for the sample choice. When relevant, describe the organism taxa, source, sex, age range and 
any manipulations. State what population the sample is meant to represent when applicable. For studies involving existing datasets, 
describe the data and its source. 


Sampling strategy Note the sampling procedure. Describe the statistical methods that were used to predetermine sample size OR if no sample-size 
calculation was performed, describe how sample sizes were chosen and provide a rationale for why these sample sizes are sufficient. 


Data collection Describe the data collection procedure, including who recorded the data and how. 
Timing and spatial scale |/ndicate the start and stop dates of data collection, noting the frequency and periodicity of sampling and providing a rationale for 
these choices. If there is a gap between collection periods, state the dates for each sample cohort. Specify the spatial scale from which 


the data are taken 


Data exclusions If no data were excluded from the analyses, state so OR if data were excluded, describe the exclusions and the rationale behind them, 
indicating whether exclusion criteria were pre-established. 


Reproducibility Describe the measures taken to verify the reproducibility of experimental findings. For each experiment, note whether any attempts to 
repeat the experiment failed OR state that all attempts to repeat the experiment were successful. 


Randomization Describe how samples/organisms/participants were allocated into groups. If allocation was not random, describe how covariates were 
controlled. If this is not relevant to your study, explain why. 


Blinding Describe the extent of blinding used during data acquisition and analysis. If blinding was not possible, describe why OR explain why 
| blinding was not relevant to your study. 


Did the study involve field work? [_] Yes [ ]No 


Field work, collection and transport 


Field conditions Describe the study conditions for field work, providing relevant parameters (e.g. temperature, rainfall). 


Location State the location of the sampling or experiment, providing relevant parameters (e.g. latitude and longitude, elevation, water depth). 
Access & import/export | Describe the efforts you have made to access habitats and to collect and import/export your samples in a responsible manner and in 
compliance with local, national and international laws, noting any permits that were obtained (give the name of the issuing authority, 


the date of issue, and any identifying information). 


Disturbance Describe any disturbance caused by the study and how it was minimized. 


Reporting for specific materials, systems and methods 


We require information from authors about some types of materials, experimental systems and methods used in many studies. Here, indicate whether each material, 
system or method listed is relevant to your study. If you are not sure if a list item applies to your research, read the appropriate section before selecting a response. 


Materials & experimental systems Methods 


n/a | Involved in the study n/a | Involved in the study 
Antibodies L] ChIP-seq 
Eukaryotic cell lines L] Flow cytometry 
Palaeontology and archaeology L] MRI-based neuroimaging 


Animals and other organisms 


Clinical data 


Dual use research of concern 
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Antibodies 


Antibodies used Describe all antibodies used in the study; as applicable, provide supplier name, catalog number, clone name, and lot number. 
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Validation Describe the validation of each primary antibody for the species and application, noting any validation statements on the 
manufacturer's website, relevant citations, antibody profiles in online databases, or data provided in the manuscript. 


Eukaryotic cell lines 


Policy information about cell lines and Sex and Gender in Research 


Cell line source(s) State the source of each cell line used and the sex of all primary cell lines and cells derived from human participants or 
vertebrate models. 


Authentication Describe the authentication procedures for each cell line used OR declare that none of the cell lines used were authenticated. 


Mycoplasma contamination Confirm that all cell lines tested negative for mycoplasma contamination OR describe the results of the testing for 
mycoplasma contamination OR declare that the cell lines were not tested for mycoplasma contamination. 


Commonly misidentified lines Name any commonly misidentified cell lines used in the study and provide a rationale for their use. 
(See ICLAC register) 


Palaeontology and Archaeology 


Specimen provenance Provide provenance information for specimens and describe permits that were obtained for the work (including the name of the 
issuing authority, the date of issue, and any identifying information). Permits should encompass collection and, where applicable, 


export. 

Specimen deposition Indicate where the specimens have been deposited to permit free access by other researchers. 

Dating methods If new dates are provided, describe how they were obtained (e.g. collection, storage, sample pretreatment and measurement), where 
they were obtained (i.e. lab name), the calibration program and the protocol for quality assurance OR state that no new dates are 
provided. 


L] Tick this box to confirm that the raw and calibrated dates are available in the paper or in Supplementary Information. 


Ethics oversight Identify the organization(s) that approved or provided guidance on the study protocol, OR state that no ethical approval or guidance 
was required and explain why not. 


Note that full information on the approval of the study protocol must also be provided in the manuscript. 


Animals and other research organisms 


Policy information about studies involving animals; ARRIVE guidelines recommended for reporting animal research, and Sex and Gender in 
Research 


Laboratory animals For laboratory animals, report species, strain and age OR state that the study did not involve laboratory animals. 


Wild animals Provide details on animals observed in or captured in the field; report species and age where possible. Describe how animals were 
caught and transported and what happened to captive animals after the study (if killed, explain why and describe method; if released, 
say where and when) OR state that the study did not involve wild animals. 


Reporting on sex Indicate if findings apply to only one sex; describe whether sex was considered in study design, methods used for assigning sex. 
Provide data disaggregated for sex where this information has been collected in the source data as appropriate; provide overall 
numbers in this Reporting Summary. Please state if this information has not been collected. Report sex-based analyses where 
performed, justify reasons for lack of sex-based analysis. 


Field-collected samples |For laboratory work with field-collected samples, describe all relevant parameters such as housing, maintenance, temperature, 
photoperiod and end-of-experiment protocol OR state that the study did not involve samples collected from the field. 


Ethics oversight Identify the organization(s) that approved or provided guidance on the study protocol, OR state that no ethical approval or guidance 
was required and explain why not. 


Note that full information on the approval of the study protocol must also be provided in the manuscript. 


Clinical data 


Policy information about clinical studies 
All manuscripts should comply with the ICMJE guidelines for publication of clinical research and a completed CONSORT checklist must be included with all submissions. 


Clinical trial registration | Provide the trial registration number from ClinicalTrials.gov or an equivalent agency. 
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Study protocol Note where the full trial protocol can be accessed OR if not available, explain why. 
Data collection Describe the settings and locales of data collection, noting the time periods of recruitment and data collection. 
Outcomes Describe how you pre-defined primary and secondary outcome measures and how you assessed these measures. 


Dual use research of concern 


Policy information about dual use research of concern 


Hazards 


Could the accidental, deliberate or reckless misuse of agents or technologies generated in the work, or the application of information presented 
in the manuscript, pose a threat to: 


No | Yes 


[ TIL] Public health 
L] E National security 
L] E Crops and/or livestock 


L] L] Ecosystems 


O [] Any other significant area 


Experiments of concern 


Does the work involve any of these experiments of concern: 


Demonstrate how to render a vaccine ineffective 


Confer resistance to therapeutically useful antibiotics or antiviral agents 


Enhance the virulence of a pathogen or render a nonpathogen virulent 


Alter the host range of a pathogen 


Enable evasion of diagnostic/detection modalities 


OO 
[TE] 
[TE] 
L] E Increase transmissibility of a pathogen 
[TE] 
[TE] 
LT] 


Enable the weaponization of a biological agent or toxin 


L] L] Any other potentially harmful combination of experiments and agents 


ChIP-seq 


Data deposition 
L] Confirm that both raw and final processed data have been deposited in a public database such as GEO. 


L] Confirm that you have deposited or provided access to graph files (e.g. BED files) for the called peaks. 


Data access links For "Initial submission" or "Revised version" documents, provide reviewer access links. For your "Final submission" document, 
May remain private before publication. | provide a link to the deposited data. 


Files in database submission Provide a list of all files available in the database submission. 


Genome browser session 
(e.g. UCSC) 


Provide a link to an anonymized genome browser session for "Initial submission" and "Revised version" documents only, to 
enable peer review. Write "no longer applicable" for "Final submission" documents. 


Methodology 


Replicates Describe the experimental replicates, specifying number, type and replicate agreement. 


Sequencing depth Describe the sequencing depth for each experiment, providing the total number of reads, uniquely mapped reads, length of reads and 
whether they were paired- or single-end. 


Antibodies Describe the antibodies used for the ChIP-seq experiments; as applicable, provide supplier name, catalog number, clone name, and lot 
number. 


Peak calling parameters | Specify the command line program and parameters used for read mapping and peak calling, including the ChIP, control and index files 
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used. 
Data quality Describe the methods used to ensure data quality in full detail, including how many peaks are at FDR 5% and above 5-fold enrichment. 
Software Describe the software used to collect and analyze the ChIP-seq data. For custom code that has been deposited into a community 


repository, provide accession details. 


Flow Cytometry 


Plots 


Confirm that: 
LI The axis labels state the marker and fluorochrome used (e.g. CD4-FITC). 


L] The axis scales are clearly visible. Include numbers along axes only for bottom left plot of group (a 'group' is an analysis of identical markers). 


L] All plots are contour plots with outliers or pseudocolor plots. 


L] A numerical value for number of cells or percentage (with statistics) is provided. 


Methodology 


Sample preparation Describe the sample preparation, detailing the biological source of the cells and any tissue processing steps used. 
Instrument Identify the instrument used for data collection, specifying make and model number. 
Software Describe the software used to collect and analyze the flow cytometry data. For custom code that has been deposited into a 


community repository, provide accession details. 


Cell population abundance Describe the abundance of the relevant cell populations within post-sort fractions, providing details on the purity of the 
samples and how it was determined. 


Gating strategy Describe the gating strategy used for all relevant experiments, specifying the preliminary FSC/SSC gates of the starting cell 
population, indicating where boundaries between "positive" and "negative" staining cell populations are defined. 


L] Tick this box to confirm that a figure exemplifying the gating strategy is provided in the Supplementary Information. 


Magnetic resonance imaging 


Experimental design 


Design type Indicate task or resting state; event-related or block design. 


Design specifications Specify the number of blocks, trials or experimental units per session and/or subject, and specify the length of each trial 
or block (if trials are blocked) and interval between trials. 


Behavioral performance measures — | State number and/or type of variables recorded (e.g. correct button press, response time) and what statistics were used 
to establish that the subjects were performing the task as expected (e.g. mean, range, and/or standard deviation across 
subjects). 


Acquisition 


Imaging type(s) Specify: functional, structural, diffusion, perfusion. 
Field strength Specify in Tesla 
Sequence & imaging parameters Specify the pulse sequence type (gradient echo, spin echo, etc.), imaging type (EPI, spiral, etc.), field of view, matrix size, 


slice thickness, orientation and TE/TR/flip angle. 


Area of acquisition State whether a whole brain scan was used OR define the area of acquisition, describing how the region was determined. 


Diffusion MRI |. ] Used |] Not used 


Preprocessing 


Preprocessing software Provide detail on software version and revision number and on specific parameters (model/functions, brain extraction, 
segmentation, smoothing kernel size, etc.). 
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Normalization If data were normalized/standardized, describe the approach(es): specify linear or non-linear and define image types used for 
transformation OR indicate that data were not normalized and explain rationale for lack of normalization. 


Normalization template Describe the template used for normalization/transformation, specifying subject space or group standardized space (e.g. 
original Talairach, MNI305, ICBM 152) OR indicate that the data were not normalized. 


Noise and artifact removal Describe your procedure(s) for artifact and structured noise removal, specifying motion parameters, tissue signals and 
physiological signals (heart rate, respiration). 


Volume censoring Define your software and/or method and criteria for volume censoring, and state the extent of such censoring. 


Statistical modeling & inference 


Model type and settings Specify type (mass univariate, multivariate, RSA, predictive, etc.) and describe essential details of the model at the first and 
second levels (e.g. fixed, random or mixed effects; drift or auto-correlation). 


Effect(s) tested Define precise effect in terms of the task or stimulus conditions instead of psychological concepts and indicate whether 
ANOVA or factorial designs were used. 


Specify type of analysis: | ] Whole brain [| ] ROl-based [ ] Both 


Statistic type for inference Specify voxel-wise or cluster-wise and report all relevant parameters for cluster-wise methods. 
(See Eklund et al. 2016) 


Correction Describe the type of correction and how it is obtained for multiple comparisons (e.g. FWE, FDR, permutation or Monte Carlo). 


Models & analysis 


n/a | Involved in the study 


L] L] Functional and/or effective connectivity 


L] L] Graph analysis 


| | Multivariate modeling or predictive analysis 


Functional and/or effective connectivity Report the measures of dependence used and the model details (e.g. Pearson correlation, partial correlation, 
mutual information). 


Graph analysis Report the dependent variable and connectivity measure, specifying weighted graph or binarized graph, 
subject- or group-level, and the global and/or node summaries used (e.g. clustering coefficient, efficiency, 
etc.). 


Multivariate modeling and predictive analysis | Specify independent variables, features extraction and dimension reduction, model, training and evaluation 
metrics. 


